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Abstract.  Despite  the  fact  that  the  main  causes  of  the  differences 
between  the  observed  Earth  nutation  and  that  derived  from  analytical 
calculations  come  from  geophysical  effects  associated  with  nonrigidity 
(core  flattening,  core-mantle  interactions,  oceans,  etc...),  efforts  have  been 
made  recently  to  compute  the  nutation  of  the  Earth  when  it  is  considered 
to  be  a  rigid  body,  giving  birth  to  several  “rigid  Earth  nutation  models.” 
The  reason  for  these  efforts  is  that  any  coefficient  of  nutation  for  a  re¬ 
alistic  Earth  (including  effects  due  to  nonrigidity)  is  calculated  starting 
from  a  coefficient  for  a  rigid-Earth  model,  using  a  frequency-dependent 
transfer  function.  Therefore  it  is  important  to  achieve  high  quality  in  the 
determination  of  rigid-Earth  nutation  coefficients,  in  order  to  isolate  the 
nonrigid  effects  still  not  well-modeled. 

After  reviewing  various  rigid-Earth  nutation  models  which  have  been 
established  recently  and  their  relative  improvement  with  respect  to  older 
ones,  we  discuss  their  specifics  and  their  degree  of  agreement. 


1.  Introduction 

The  precision  in  the  determination  of  the  Earth’s  nutation  through  up-to-date 
techniques  such  as  Very  Long  Baseline  Interferometry  (VLBI)  has  drastically  im¬ 
proved  in  the  recent  years,  so  that  it  has  become  usual  now  to  express  corrections 
of  the  nutation  coefficients  in  terms  of  microarcseconds,  instead  of  milliarcsec- 
onds,  as  was  the  case  at  the  eve  of  the  90s.  Comparing  the  observational  data  of 
nutation  at  this  new  level  of  precision,  with  the  theoretical  determinations,  is  of 
primary  interest  for  the  knowledge  of  our  planet.  Indeed  the  way  by  which  the 
Earth  is  reacting  to  the  action  of  the  other  celestial  bodies  (the  Moon,  the  Sun 
and  the  planets)  and  thus  by  which  its  nutation  is  modified  by  external  forcing 
gives  precious  information  on  the  mechanisms  or  parameters  which  cannot  be 
so  well  known  through  other  studies. 

The  present  conventional  nutation  series  adopted  as  the  IAU  1980  nutation 
series  (Seidelmann,  1982),  although  being  the  standard  of  reference,  has  been 
found  to  differ  significantly  from  observational  data  obtained  by  VLBI  analyses. 
For  many  nutation  coefficients,  the  difference  between  the  estimated  value  from 
VLBI  analysis  and  the  analytical  one  is  much  larger  than  the  uncertainties  of 
the  data.  Following  that  fact,  an  empirical  series  of  nutation  for  practical  use 
giving  close  agreement  to  the  data  has  been  elaborated  and  is  now  available  in 
the  IERS  Conventions  (McCarthy,  1996). 
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Currently  active  research  is  devoted  to  the  interpretation  of  the  differences 
between  the  observed  nutation  and  its  amount  as  given  by  the  series,  and  to  the 
construction  of  series  starting  from  a  model  as  representative  as  possible,  of  the 
various  geophysical  effects  on  the  nutations.  The  basic  model  of  nutation  of  Wahr 
(1981),  applied  to  an  oceanless,  elastic,  ellipsoidal  Earth  model  derived  from  the 
Earth  model  1066A  of  Gilbert  and  Dziewonski  (1975)  under  the  assumption  of 
hydrostatic  equilibrium  is  no  longer  adequate.  An  accurate  model  of  nutation 
must  include  various  effects  such  as  mantle  anelasticity,  effects  of  ocean  tides  and 
electromagnetic  couplings  produced  at  the  core-mantle  boundary  (CMB)  or  the 
inner-core  boundary  (ICB),  as  was  shown  by  Mathews  et  al.  (1991).  Moreover, 
non-hydrostatic  equilibrium  has  to  be  considered  (Dehant  and  Defraigne,  1997; 
Schastok,  1997). 

In  parallel,  the  elaboration  of  these  models  requires  high  quality  nutation 
coefficients  of  the  Earth  considered  as  a  rigid  body,  which  is  a  more  straightfor¬ 
ward  and  simple  approximation.  In  that  frame  the  challenge  is  quite  clear.  It  is 
to  model  what  should  be  the  motion  of  nutation  (and  precession)  of  the  Earth 
considered  as  a  rigid  body  undergoing  the  perturbing  effects  of  the  Moon,  the 
Sun  and  the  planets,  all  considered  as  point  mass  bodies.  Their  distances  are  a 
priori  too  large  for  their  forms  to  have  significant  effects  on  the  motion.  This 
is  a  necessary  primary  step  when  one  wants  to  compute  the  coefficients  of  the 
nutation  of  the  Earth  taking  into  account  various  effects  related  to  the  oceans, 
the  atmosphere,  and  all  the  geophysical  influences  which  have  been  pointed  out 
in  the  last  paragraph.  Thus,  until  now,  all  the  theories  giving  birth  to  the 
computation  of  the  coefficients  of  nutation  (Wahr,  1981;  Dehant  and  Defraigne, 
1997;  Mathews  et  al. ,  2000;  Schastok,  1997;  Getino  and  Ferrandiz,  1999)  start 
from  the  value  of  these  same  coefficients  where  the  Earth  is  considered  as  a 
rigid  body.  Generally  the  way  to  convert  from  a  nutation  coefficient  for  a  rigid- 
Earth  model,  to  the  same  coefficient  for  a  nonrigid  Earth,  is  through  a  transfer 
function,  generally  frequency-dependent  that  includes  various  parameters  (Love 
numbers,  ellipticity  of  the  core,  dynamical  ellipticity  of  the  Earth)  mainly  of 
geophysical  nature. 

Thus  we  understand  easily  why  to  determine  the  coefficients  of  nutation  for 
a  rigid-Earth  model  is  fundamental  to  reveal  in  the  best  way  the  contributions 
arising  from  the  nonrigid  part  associated  with  the  Earth’s  interior.  Moreover, 
the  remarkable  improvement  of  the  accuracy  of  VLBI  observations  from  the  be¬ 
ginning  of  the  80s  has  led  to  an  increasing  number  of  studies  devoted  to  the 
determination  of  rigid-Earth  nutation.  These  studies  were  in  particular  moti¬ 
vated  by  the  fact  that  the  truncation  of  the  coefficients  of  nutation  of  Kinoshita’s 
series  for  a  rigid  Earth  model  (Kinoshita,  1977)  that  serves  as  the  basic  series  for 
the  IAU1980  series  of  nutation  (Seidelmann,  1982),  was  0.1  mas.  This  leads  to  a 
lack  of  precision  at  the  level  of  1  milliarcsecond  for  the  total  nutation  amplitude. 
This  shows  how  it  became  necessary  to  push  further  the  level  of  truncation  of 
the  coefficients  of  rigid-Earth  nutation. 

Presently  the  sub-microarcsecond  level  has  been  adopted  by  authors  dealing 
with  the  topic  (Roosbeek  and  Dehant,  1997;  Bretagnon  et  al .,  1999;  Souchay  et 
al .,  2000).  A  convergence  between  the  rigid-Earth  nutation  as  given  by  these 
authors  should  lead  to  the  conclusion  that  at  least  we  can  be  comforted  with  the 
validity  of  such  a  basic  part.  In  the  following  we  show  that  this  convergence  is 
really  satisfactory,  after  summarizing  the  evolution  of  the  precision  and  accuracy 
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of  the  rigid-Earth  nutation  theory  through  the  various  works  which  have  been 
devoted  to  the  matter  in  the  last  decade. 


2.  Evolution  of  the  rigid-Earth  nutation  theory  in  the  90s 

Kinoshita  and  Souchay  (1990)  published  new  series  of  nutation  for  a  rigid-Earth 
model,  considering  a  new  level  of  truncation  for  the  coefficients  of  nutation:  0.005 
mas.  This  new  level  made  necessary  the  study  of  new  effects  and  the  extension 
of  some  developments.  The  principal  influences  can  be  briefly  summed  up  as 
(Kinoshita  &  Souchay,  1990): 

•  Influence  of  the  J2  harmonic  and  of  the  second-order  parts  of  the  potential 
of  the  Earth  (J3,  triaxiality  and  J4), 

•  Direct  and  indirect  planetary  effects, 

•  Extension  of  the  theory  to  second  order:  coupling  effect  between  the  rota¬ 
tional  motion  of  the  Earth  and  the  orbital  motion  of  the  Moon. 


3.  New  developments  in  rigid-Earth  nutation  theory 

The  updated  value  of  the  general  precession  in  longitude  led  to  a  change  of  the 
value  of  the  dynamical  ellipticity  of  the  Earth  as  well  as  the  scaling  factors  for 
the  nutation  coefficients  relative  to  the  potential  of  the  Moon  and  the  Sun.  For 
this  reason,  Souchay  and  Kinoshita  (1996,  1997)  recalculated  the  coefficients 
of  nutation  —  due  to  the  lunisolar  inffuence  on  the  geopotential  (J2,  J3,  C22, 
S22  an(i  J4)  and  those  coming  from  the  direct  torque  exerted  by  the  planets 
—  influenced  by  these  changes.  Moreover,  they  presented  new  contributions  to 
the  nutation  not  included  in  Kinoshita  &  Souchay  (1990)  such  as  the  influence 
of  the  periodic  oscillations  of  the  ecliptic,  and  the  planetary  tilt-effect.  The 
agreement  of  the  new  coefficients  of  Souchay  &  Kinoshita  (1996,  1997)  with 
those  of  Hartmann  &  Soffel  (1994)  and  Williams  (1994,  1995)  was  remarkable. 


4.  The  effects  on  nutation  of  the  non-zonal  harmonics  of  third  and 
fourth  degree 

Taking  into  account  the  new  accuracy  of  the  precession-nutation  observations, 
that  is  to  say  a  few  microar seconds,  Folgueira  et  al.  (1998a, b)  calculated  the 
coefficients  of  the  nutation  for  a  rigid-Earth  model  due  to  C^m  and 

Sim  ( to  7^  0)  harmonics  of  the  geopotential,  which  were  not  considered  in  previ¬ 
ous  studies  of  nutation  involving  Hamiltonian  theory,  with  a  level  of  truncation 
of  0.1  /4as  (microarcseconds). 

These  non-zonal  harmonics,  together  with  C 22  and  622,  give  birth  to  the 
short-period  nutations,  the  periods  being  related  to  the  order  m  of  the  corres¬ 
ponding  harmonic,  that  is  to  say: 

•  the  diurnal  terms  come  from  C31,  631,  C41  and  ^i, 
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•  the  semidiurnal  terms  are  due  to  the  harmonics  with  m  =  2,  and 

•  the  terdiunal  terms  have  their  origin  in  C 33,  S 33,  C43  and  N43  harmonics. 

The  principal  diurnal  and  subdiurnal  terms  for  the  nutation  in  longitude 
and  obliquity  are  shown  in  Tables  1  and  2. 


Table  1.  Principal  terms  for  quasidiurnal  nutations  in  longitude  and 
obliquity  for  the  figure  axis.  The  unit  is  gas. 


Argument  Period  Longitude  (A ip)  Obliquity  (Ae) 


<p 

L  M 

L> 

sin 

COS 

sin 

COS 

r 

0 

0 

r 

0 

1 

0.962r5 

-38.23r3 

-4.6980 

-4.8567 

45.4063 

r 

0 

0 

-r 

0 

-r 

r.03505 

-35.4036 

-4.3543 

-4.5868 

42.9409 

r 

-r 

0 

r 

0 

r 

0.99696 

24.5443 

2.9675 

4.4636 

-9.4678 

r 

r 

0 

-r 

0 

-r 

0.99758 

-49.9400 

-2.4507 

-0.9746 

7.9055 

r 

r 

0 

r 

-2 

r 

0.992r6 

-7.0629 

-0.8684 

-0.3452 

2.8089 

Table  2.  Principal  terms  for  subdiurnal  nutations  in  longitude  and 
obliquity  for  the  figure  axis.  The  unit  is  gas. 


Argument  Period 

Longitude  (A ip) 

Obliquity  (Ae) 

O'  Im  P  D  ^ 

sin  cos 

sin  cos 

2 

0 

0 

0 

0 

0 

0.49863 

34.2522 

-47.9397 

-7.4457 

-42.3964 

2 

0 

0 

-2 

0 

-2 

0.54753 

-25.5274 

44.6533 

5.6775 

9.8906 

2 

0 

0 

-2 

2 

-2 

0.50000 

-40.4455 

5.9960 

2.3682 

4.4255 

2 

-r 

0 

-2 

0 

-2 

0.52743 

-5.0979 

2.9264 

4.4848 

2.0588 

2 

0 

0 

-2 

0 

-r 

0.54756 

-4.7405 

2.7040 

4.0846 

4.8842 

5.  The  sub-microarcsecond  rigid-Earth  nutation  series  REN  2000 

Souchay  et  al.  (1999)  presented  the  new  tables  REN  2000  of  the  nutation  of  a 
rigid  Earth,  starting  from  Hamiltonian  theory,  with  a  sub-microarcsecond  level 
of  truncation.  All  the  previously  studied  influences,  together  with  the  second- 
order  effects  due  to  crossed-nutations  and  spin-orbit  coupling,  are  included  in 
this  paper. 


6.  The  diurnal  and  sub-diurnal  nutations  of  REN  2000 

Recently,  other  rigid-Earth  nutation  series,  calculated  using  different  analytical 
methods,  have  come  forth.  These  are  SMART97  and  RDAN97  series  (Bretagnon 
et  al.,  1997,  1998;  Roosbeek  &  Dehant,  1998).  Folgueira  &  Souchay  (1999)  and 
Folgueira  et  al.  (1999)  have  compared  the  nutations  given  by  REN  2000  series 
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with  the  corresponding  ones  in  SMART97  and  RDAN97.  Table  3  gives  the 
comparison  for  the  short-period  nutations  in  the  time  domain. 

Moreover,  they  have  computed  corresponding  values  for  a  nonrigid  Earth 
of  REN  2000  prograde  nutations,  using  a  new  transfer  function  developed  by 
Mathews  (1999)  and  showed  that  some  coefficients  changed  at  the  level  of  a  few 
lias.  Finally,  they  have  also  explained  how  short-period  nutations  can  perturb 
the  analysis  of  the  polar  motion  (Bizouard  et  al. ,  1999a;b). 


Table  3.  Comparison  in  time  domain  between  SMART97,  RDAN97 
and  REN  2000  short  period  nutations.  The  unit  is  jias. 


Aip  sin  So 

SMART97  -  REN  2000 

1.74 

3.52 

RDAN97  -  REN  2000 

2.40 

6.80 

RDAN97  -  SMART97 

1.40 

6.80 

7.  Conclusions 

In  the  present  paper,  we  have  analyzed  the  different  steps  carried  out,  in  this 
last  decade,  to  reach  at  the  final  sub-microarcsecond  tables  of  nutation  for  a 
rigid-Earth  model  REN  2000  (Souchay  et  al. ,  1999).  We  can  conclude  that: 

1.  The  tables  REN  2000  of  the  nutation  for  a  rigid  Earth  model,  starting  from 
the  Hamiltonian  theory,  catch  all  the  coefficients  at  a  sub-microarcsecond 
level. 

2.  There  is  a  good  agreement  between  REN  2000  and  other  recent  works 
about  rigid-Earth  nutations  (Bretagnon  et  al .,  1997,  1998  ;  Roosbeek  & 
Dehant,  1998). 

3.  The  short-period  nutations  are  the  last  most  important  contribution  to 
reach  at  the  0.1  /ias  level  in  the  REN  2000  tables. 

4.  The  validity  of  the  analytical  series  REN  2000  has  been  checked  by  a 
numerical  integration  of  the  nutation  developed  by  Souchay  &  Kinoshita 
(1991)  and  Souchay  (1998)  and  by  the  study  of  the  residuals  between  the 
results  given  by  these  two  different  methods. 
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